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ABSTRACT 


A  ray  theory  for  the  transmission  of  light  across  a  ran¬ 
domly  rough  boundary  is  formulated.  Multiple-scatter  and 
shadowing  effects  are  accounted  for  in  a  consistent  way. 
The  illumination  probabilities  for  the  transmission  prob¬ 
lem  are  derived  in  a  nontrivial  extension  of  current  shad 
owing  theory.  The  direct  transmission  term  and  the  two 
main  higher-order  corrections  are  specified  in  complete 
detail  for  a  normally  distributed  boundary.  The  ap¬ 
plication  of  transmission  theory  to  a  quantitative  des¬ 
cription  of  sea  color  is  considered. 
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I .  INTRODUCTION 

A  comprehensive  study  of  t.ic  sea  requires  investigations  in  both  the 
microwave  and  optical  portions  of  the  spectrum.  Thus,  microwave  radiometry 
provides  a  sensitive  all-weather  tool  for  the  calculation  of  boundary  tempera¬ 
tures  and  sea  states,  while  subsurface  properties  can  be  probed  by  the  spectral 
variation  of  sea  color.  For  example,  the  chlorophyll  content  in  the  ocean, ^ 
the  depth  of  coastal  waters,  and  the  degree  of  pollution  in  coastal  and  in¬ 
land  waters  all  may  be  identifiable  by  their  spectral  signatures. 

The  problem  of  sea  color  is  very  complex.  Sunlight  and  skylight  are  both 
scattered  by  the  sea  surface  and  transmitted  across  the  boundary  into  the  water 
medium.  Some  of  this  transmitted  energy  is  backscattered,  either  by  molecular 
and  particulate  matter  in  deep  water  or  by  the  bottom  in  shallow  waters,  so 
that  a  certain  fraction,  the  upwelling  sea  light,  is  transmitted  back  across 
the  boundary  into  the  air.  Sea-color  measurements  must  be  made  at  angles  out¬ 
side  the  glitter  pattern  of  the  scattered  sunlight,  so  the  sea  color  is  deter¬ 
mined  by  the  two  components  of  scattered  skylight  and  upwelling  sea  light. 

The  actual  analysis  of  the  sea-color  problem  must  account  for  surface 
roughness.  Thus,  sunlight  incident  on  a  flat  ocean  surface  will  be  transmitted 
in  the  single  direction  determined  by  Snell's  law.  For  the  general  ocean 
boundary,  however,  Snell's  law  must  be  applied  to  each  of  the  locally  flat 
elements  compris'ng  the  surface.  This  results  in  an  angular  distribution  of 
transmitted  energy ,  which,  of  course,  will  produce  a  final  component  of  up- 
welling  sea  light  with  angular  distribi. jion  dependent  on  the  surface  sea  state. 

The  scattered  skylight  will  also  have  an  angular  radiation  pattern  peculiar  to 
the  sea  state.  Quite  clearly,  if  subsurface  properties  are  of  primary  interest, 
we  want  to  make  sea-color  observations  in  directions  where  the  scattered  sky¬ 
light  does  not  dominate  the  upwelling  sea  light.  Thus,  experiments  in  this 
field  cannot  be  designed,  and  the  data  properly  interpreted,  without  an  analy¬ 
tical  description  of  the  interaction  of  optical  radiation  with  irregular  surfaces. 

Until  recently,  the  theory  of  scattering  from  rough  surfaces  was  severly 
limited,  even  in  the  geometrical -optics  realm  appropriate  to  optical  studies. 

This  was  due  to  an  incomplete  treatment  of  the  geometri cal -opti cs  description, 
in  which  the  nonlocal  effects  of  shadowing  and  multiple  scatter  of  the  rays 
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were  neglected.  This  neglect  was  manifested  in  the  inability  of  the  theory  to 
conserve  energy.^  An  extended  scattering  theory  now  exists^  which  accounts  for 
shadowing  and  double-scatter  effects  in  a  fully  consistent  way,  i.e.,  energy  is 
conserved  to  a  high  degree  of  accuracy.  While  all  the  tools  necessary  for  an 
accurate  determination  of  the  scattered  skylight  are  now  at  hand,  there  has 
been  no  corresponding  effort  directed  towards  a  theory  of  rough-surface  trans¬ 
mission.  In  this  report,  we  present  a  complete  geometrical -optics  theory  of 
transmission,  with  shadowing  and  multi pie- reflection  effects  entering  the  des¬ 
cription  in  a  natural  and  consistent  way.  There  are  no  restrictions  on  incidence 
or  observation  angles,  and  the  full  range  of  ocean  sea  states  can  be  treated 
accurately. 

We  treat  the  idealized  problem  of  two  different  semi -infinite  and  homo¬ 
geneous  media  separated  by  a  normally  distributed  rough  boundary.  A  source  is 
located  in  one  medium,  and  the  observer  is  located  in  the  other.  In  Section  II, 
Vie  define  the  theoretical  quantities  of  interest,  namely  the  transmission  co¬ 
efficient  which  describes  the  angular  radiation  pattern  after  refraction,  and 
the  absorptivity,  which  gives  the  total  fraction  of  incident  energy  entering 
the  observer  medium.  In  Section  III,  a  formal  analysis  of  the  absorptivity  is 
presented  in  terms  cf  the  physical  subprocesses  into  which  energy  may  be  chan¬ 
neled,  Thus,  the  dominant  process  would  be  simple  absorption  of  an  incident 
ray,  but,  for  a  rough  boundary,  we  can  also  have  absorption  followed  by  re¬ 
flection  bivlow  the  surface,  as  well  as  multiple  reflection  above  the  surface 
followed  by  absorption.  The  description  is  further  complicated  by  the  small, 
but  nonzero,  probability  that  a  light  ray  may  be  transmitted  into  an  undulation 
of  the  observer  medium  only  to  be  refracted  back  into  the  source  medium.  A 
consistent  approach  to  the  cumulative  effects  of  all  the  various  higher-order 
Processes  is  presented. 

Section  IV  contains  an  analogous,  and  very  formal , representation  of  the 
second  theoretical  quantity,  the  transmission  coefficient,  in  terms  of  physical 
subprocesses.  The  first  actual  application  of  the  transmission  theory  to 
random  rough  surfaces  is  carried  through  in  Section  V.  Specifically,  we  derive 
the  direct  transmission  coefficient  which  describes  the  angular  transmitted 
radiation  pattern  in  the  absence  of  multiple-reflection  contributions.  Shadowing 
corrections  are  incorporated  automatically.  A  simple  exponential  model  is  used 
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to  account  for  propagation  losses.  The  result  is  applicable  to  the  full  two 
dimensional  rough  surface  and,  as  such,  it  contains  all  polarization  informa¬ 
tion.  Not  surprisingly,  the  strength  of  the  transmission  at  a  given  observation 
angle  is  proportional  to  the  probability  density  of  that  slope  which  connects 
the  incident  and  observation  directions  by  Snell's  law.  It  is  of  interest  that 
the  direct  transmission  coefficient  is  identically  zero  for  some  ranges  of  ob¬ 
servation  angle,  regardless  of  the  rms  surface  slope. 

The  results  of  a  shadowing  theory  tailored  to  the  transmission  problem 
are  merely  cited  in  Section  V,  but  in  Section  VI  we  present  a  full  analytical 
treatment.  The  illumination  probabilities  derived  here  are  not  trivial  exten¬ 
sions  of  the  illumination  probabilities  used  in  scattering  theory.  This  follows 
because  the  source  and  observer  are  now  on  opposite  sides  of  the  rough  boundary, 
so  that  a  high  point  on  the  surface  relative  to  the  source  is  a  low  point  on 
the  surface  relative  to  the  observer.  The  results  are  obtained  in  terms  of 
gamma  functions,  which  are  fully  tabulated. 

higher-order  transmission  terms  are  formulated  in  Sections  VII  and  VIII 
for  the  mathematically  tractable  one  dimensional  rough  surface  (cylindrical 
symmetry).  Thus,  a  ray  of  sunlight  may  reflect  from  one  point  on  the  ocean 
surface  only  to  Intersect  the  surface  again  at  another  point.  There  will  be 
transmission  across  the  surface  boundary  at  the  second  point  as  well  as  direct 
transmission  at  the  first  scatter  point.  The  angular  profile  due  to  trans¬ 
mission  at  the  second  point  fills  all  possible  observation  angles,  but  the  total 
fraction  of  incident  energy  associated  with  this  double-scatter  process  is 
estimated  at  only  about  3%  of  the  directly  transmitted  energy  for  likely  ocean 
roughnesses.  The  corresponding  process,  in  which  a  directly  transmitted  ray 
of  sunlight  is  internally  reflected  by  the  rough  boundary  is  also  formulated. 

It  is  expected  to  contribute  significantly  to  the  below-surface  radiation 
profile  only  for  near-grazing  incidence  and  the  very  roughest  of  ocean  surfaces. 

In  Section  IX,  we  consider  the  inverse  problem  of  refraction  into  the 
medium  of  lower  index  (e.g.,  water-to-air).  Fortunately,  it  is  not  necessary 
to  repeat  the  previous  calculations,  but  instead  we  apply  the  principle  of 
reciprocity.  This  theorem  relates  the  effects  of  an  interchange  of  source  and 
observer,  and  so  is  ideally  suited  for  application  to  the  inverse  problem. 


3 


Finally,  in  Section  X,  we  derive  the  formal  relation  between  observed  intensity 
and  the  transmission  and  scattering  coefficients  which  must  be  used  in  any 
analytical  treatment  of  sea  color.  A  homogeneous  sea  with  flat  bottom  (of 
zero  slope)  is  considered  as  a  special  case. 

II.  PROBLEM  STATEMENT  AND  DEFINITIONS 

We  consider  the  idealized  situation  of  two  differ-ent  homogeneous  and 
semi-infinite  media  separated  by  an  irregular  boundary  (Fig.  1).  One  of  the 
media  contains  a  source,  with  the  receiver  located  in  the  second  medium.  For 
example,  sunlight  is  incident  on  a  rough  ocean  surface  and  we  are  interested 
in  tne  angular  distribution  of  transmitted  radiation  below  the  surface.  We 
will  take  \:he  index  of  refraction  of  the  source  medium  as  unity,  whereas  the 
second  medium  will  be  represented  everywhere  by  a  complex  constant  n.  The 
homogeneity  assumption  means  that  once  clear  of  the  surface,  the  refracted 
rays  propagate  recti  linearly  until  completely  absorbed.  A  simple  exponential 
absorption  is  associated  with  each  ray.  We  will  also  restrict  the  investi- 
gation  to  those  cases  where  the  imaginary  part  of  the  refractive  index  is 
small,  e.g.,  the  problem  of  optical  transmission  i.ito  clear  water.  Then 
Snell's  law  relates,  with  high  accuracy,  the  angle  of  ref» action  to  n,  the 
real  part  of  n  (Appendix  A). 

There  are  two  theoretical  quantities  of  interest  to  us,  The  first  is 
the  absorptivity,  defined  as  the  fraction  of  incident  energy  which  is  not 
scattered.  The  absorptivity  is  a  function  only  of  the  angle  of  incidence. 

Now,  energy  crossing  the  boundary  must  propagate  some  distance  into  the  medium 
(of  index  n)  before  it  is  absorbed.  This  view  leads  immediately  to  the  second 
theoretical  quantity,  the  transmission  coefficient,  which  determines  the  frac¬ 
tion  of  incident  energy  observed  at  some  angle  relative  to  the  (-z)  axis  and 
at  a  depth  Z,  Thus,  the  transmission  coefficient  describes  the  angular  spread¬ 
ing,  due  to  surface  roughness,  of  the  refracted  radiation,  and  it  also  accounts 
for  propagation  losses.  If  we  designate  a(9Q»<l>Q)  at  the  absorptivity  due  to  a 
source  at  (6^,(1)^),  and  if  t(e,(j);  is  the  transmission  coefficient 

for  an  observer  at  (0,(())  and  at  depth  £  in  a  medium  with  absorption  coefficient 
K,  then 

a(3o><t>o)  =  ®o’'^o* 

whera  the  integration  is  over  the  lower  hemisphere  (all  possible  refracted  angles). 
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III.  PARTIAL-SURFACE  REPRESENTAFION  FOR  ABSORPTIVITY 

Radiation  incident  on  a  perfectly  flat  surface  is  refracted  according  to 
Snell's  law.  For  a  rough  ocean  surface  and  optical  radiation,  we  can  apply  Snell's 
law  to  each  of  the  locally  flat  surface  elements  comprising  the  boundary.  Optical 
wavelengths  will  be  small  compared  with  even  capillary  structure  on  the  ocean  sur¬ 
face,  so  a  ray  theory  of  transmission  can  be  expected  to  give  very  accurate  results. 

A  ray  theory  of  transmission  for  an  irregular  surface  is  actually  the  sum  of 
many  processes.  The  most  important  are  illustrated  by  Figs.  2-4.  Every  incident  ray 
contributes  to  the  direct  transmission  process,  (Fig»  2).  Figure  3  illus¬ 
trates  an  additional  transmission  contribution,  reflec¬ 

tion  of  the  incident  ray  at  point  1  and  a  subsequent  intersection  with  the  surface 
at  point  2.  The  importance  of  this  process  liicrcases  with  increasing  rms  surface 
slope.  The  associated  process  of  a  transmitted  ray  which  undergoes  a  subsequent 
reflection  (below  the  surface)  at  point  2,  kj^kjy>-t2R»  ‘'s  shown  in  Fig.  4.  A  cor¬ 
rect  theory  of  transmission  must  account  for  all  such  diagrams. 

In  this  section  we  will  restrict  the  quantitative  discussion  to  the  special 
one-d’mensional  rough  surface  (cylindrical  symmetry)  in  order  to  emphasize  the 
physical  picture;  cross  polarization  complications  are  introduced  in  Section  V. 

Thus,  we  consider  radiation  of  unit  intensity  and  beam  width  incident  on  an  ir¬ 
regular  surface  described  by  z  =  c(x,y)  =  ?(x).  The  incident,  transmitted,  and 
scattered  rays  are  all  contained  in  the  x=z  plane.  The  optical  wavelength  is 
small  compared  to  all  surface  parameters.  The  scattering  and  transmission  pro¬ 
cesses  are  incoherent  in  this  geometrical -optics  realm  as  each  surface  element 
simply  obeys  Snell's  law,  independent  of  any  other  surface  element. 

We  imagine  a  ray  trace  for  each  surface  element,  and  we  will  include 
that  element  in  one  of  a  number  of  classes,  depending  on  the  scattering  char¬ 
acter  of  that  trace.  Specifically,  we  divide  the  total  surface  into  two  parts, 

2  =  JHq',  where  every  surface  element  in  n  is  visible  to  the  incident  beam 
with  direction  k^.  The  elements  in  are  shielded  from  the  incident  beam, 
and  as  their  contributions  are  multiplied  by  zero  incident  intensity,  their  ray 
traces  have  zero  weight.  The  illuminated  subset  is  now  split  into  two  parts, 
n  =  nj+Jij,  where  every  element  in  is  characterized  by  a  specularly  scattered 
ray  which  does  not  intersect  the  surface  elsewhere  (Fig.  2).  For  every 
element  of  the  subset  ^1  intersect  the  surface 
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again*  but  now  is  defined  as  a  smaller  subset  with  the  property  that  the 
follow-up  ray  does  not  intersect  the  surface  elsewhere  (Fig.  3).  Me  con¬ 

tinue  this  bookkeeping  until  we  achieve  an  empty  subset,  say  (a  maximum 
of  m  surface  scatterings):  , 

,  s!  j  n,.  .  ,  ,  (3.1) 

The  incident  power  intercepted  by  a  surface  element  dSj^  in  n  is  unit  in¬ 
tensity  times  the' elemental  area  projected  onto  the  indident  wave  front,  i.e. , 
(-ICq* njJdSi.  The  initial  polarization  is  taken  as  either  vertical  or  hori¬ 
zontal ,  and  the  scattered  or  transmitted  radiation  will  retain  this  polarization 
because  of  the  cylindrical  symmetry.  If  the  reflectivity  at  surface  point  i , 
is  r^  (for  either  vertical  or  horizontal  polarization),  then  (1-r..)  is  the 
fraction  of  energy  incident  at  this  po1|fit  which  is  transmitted.  We  are  now  in 
a  position  to  write  down  the  quantity  '  T(ejj),  defined  as  the  total  fraction  of 
energy-  incident  from  k^  which, is  rdf racted  across  the  boundary: 


=  s. 


m 

1  -1  I 


=  |dSj(cosa2)[i-rj(cosa2j] 

Ut 


dS  j  ( cosa  j )  r  j  ( cosa  j )  1 1-  rg  ( cosctg ) ) 


■  .  ..  [dSj(cos(Sj^)rjr2 

fi  ’  1  ,  , 
m-1  ' 


’'m  i(lri‘  ) 

m-1'  m' 


(3.2) 


where  (“^q’^i)  =  cosa^  (-Cj^rig)  ?  cosa2,  and  an  obvious  abbreviation  for  the 
arguments  Qf  the  reflectivities  is  introduced  In  the  m'th  term.  Eqich  inte¬ 
gration  is  taken  over  the  indicated  subset  of  the  total  surface.  Equation 
(3.2)  accounts  for  all  orders  of  multiple  reflection  in  the  source  medium. 

For  example,  the  contribution  to  T(0jj)  at  point  1  in  Fig.  3  is  contained  in 
the  first  term,  while  the  energy  crossing  the  boundary  at  point  2  of  the  same 

'  I 

diagram  is  accounted  for  in  the  second  term, 
normal izati, on  factor,  for  when  r  =  0  and  T  =  1,  we  must  have 


Note,  the  quantity  is  the 


(3.3) 


n 


II 


/ 


.i. 


# 


I 

I 

( 

,  The  absorptivity  is  defined  as  the  fraction  of  incident  energy  which  is  not 

scattered.  Unfortunately,  the  quantity  1(9^)  is  not  the  absorptivity.  This 
follows  because  energy  refracted  across  a  rough-surface  interface  can  be  re¬ 
fracted  back' into  the  source  medium,  there  to  scatter  and  refract  again.  Figure 
4  contains  lan  illustration  of  this  point.  Thus,  an  incident  ray  refracts  at 
point  1,  then  propagates  internally  to  point  2  where  it  reflects  and  may  refract 
i(depending  on  the  slope  at  point  2)  back  into  the  source  medium.  The  doubly 
refracted  ray  kgjj  is  always  directed  towards  the  surface,  so  it  will  be  re¬ 
flected  and  refracted  again  at  point  3.  The  double  transmission  plus  reflection 
process  augments  the  scattered  radiation  pattern.  Therefore,  T  is  greater  than 
the  absorptivity  by  an  additional  scatter  contribution  6(0^): 

'  ,  T(^o)  =  aCe^j)  +  6(0^^)  .  (3.4) 

In  ordeh  to  obtain  a  partial -surface  expansion  for  3(6^^)  alone,  we  must  expand 
each  of  the  T^*  in  Eq.  (3.2)*  i  a  partial-surface  representation.  We  will  simply 
state  in  passing  that,  while  TCe^)  is  not  strictly  the  fundamental  quantity  of 
interest,  it  may  be  bounded  in  value  from  both  above  and  below.  As  6{6q)  is 
expected  to  be  very  small  in  magnitude,  this  can  be  an  important  advantage  in 
numerical  computations  of  the  absorptivity. 

I 

We  now  expand  Tj  in  a  partial -surface  series  which  accounts  for  all  the 
possible  processes  the  ray  may  undergo  after  refraction  at  point  1.  We  have, 

I  expl i ci tly , 

SqTiOo)  =  |dS^  cosa^d-rj)  ,  (3.5) 

Q 

and  we  define 

n  =  fly  , 

flj  -  ^TR  ’  (^•®) 

Here,  fly  contains  those  surface  elements  in  such  that  the  refracted  ray  kjy 
does  not  intersect  the  surface  again  (Fig.  2),  whereas  the  ray  k^y  does  inter¬ 
sect  the  surface  again  for  elements  of  surface  in  the  subset  (Fig*  4). 

I  *  * 

The  further  decomposition  of  fly  depends  on  whether  the  ray  reflected  from 
point  2  in  Fig.  4  clears  the  surface  (ily[^)  or  intercepts  it  again  (Sl-j-p). 

( 

f 

7 


/ 


I 


We  can  write,  with  these  definitions, 

SoTi(0o)  =  |dSj(cosaj)sl-rj)  +  |d$j(cosaj)(l-rj)[r2+(l-r2)] 

fly  ^ 


Fig.  4;  some  of  the  energy  associated  with  this  term  goes  into  absorption  and 
some  into  scattering. 

A  most  important  aspect  to  the  expansion  of  Eq.  (3.5)  is  that  the  value 
of  the  residue  term,  ^^(0^),  can  be  calculated  indirectly.  This  is  true  because 
Ti(eo)  is  nothing  more  than  the  direct  emissivity^  of  the  rough  surface.  Thus, 
a  computer  program  already  exists  for  the  evaluation  of  Tj(9p).  Some  results 
are  listed  in  Tables  I  and  II  for  a  rough  water  surface  and  incident  radiation 
of  wavelength  0.7vi.  In  addition,  both  a^  and  a^j^  can  be  calculated  numerically 
(although  we  have  not  yet  done  this),  so  we  can  solve  for  Aj  by  use  of  Eq.  (3.7). 
Alternatively,  we  can  limit  numerical  computations  to  T,  and  a^  and  solve  for 
the  value  of  (a^j^+Aj).  This  will  ordinarily  be  sufficient  as  a^j^  is  surely 
negligible  compared  to  a^  save  for  near  grazing  incidence  and  very  rough  sur¬ 
faces  (Section  VIII). 


where 

(■3.7) 

J()Sj(cosc«[)(l-rj)  , 

(3,8) 

jdSj(cosaj)(l-rj)r2  , 

«TR 

(3.9) 

jdSj(cosaj)(l-rj)r2  +  |dSjCosaj(l-rj)(l-r2) 

n  Q' 

TR  T 

(3.10) 

The  quantities  a^  and  a^j^  are  true  contributions  to  the  absorptivity,  for  Jlj 
and  were  constructed  such  that  and  k2pj,  respectively,  do  not  intercept 

the  surface  again.  On  the  other  hand,  Aj  contains  higher-order  absorption 
processes  as  well  as  contributions  to  the  scattering  diagram.  For  example, 
the  second  term  in  Eq.  (3.10)  is  the  double  refraction  process  illustrated  in 

m 
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Table  I 


The  quantities  T-j(6q)  and  T2(eQ),  evaluated  for  a  rough 
water  surface  of  15°  rms  slope  and  horizontally  polarized  incii 
radiation  of  wavelength  0.7y. 

®0  t^l(®0^^  ""  ^13(60)]  X  293  T2/T^  X  100 

30° 

283.0  0.47 

0.17 

50 

271.4  3.03 

1,11 

65 

252.4  7.32 

2.90 

75 

235.0  7.24 

3.08 

80 

225.5  7.36 

3.26 

85 

215,4  4.11 

Table  II 

1.91 

The  quantities  T^(6q)  and  T2{ep),  evaluated  as  in  Table  I, 
but  for  vertically  polarized  incident  radiation. 

!o 

[T^(0^,)]  X  293  [12(60)]  X  293 

T2/T^  X  100 

30° 

288.3  0.12 

0.04 

50 

287.6  1.53 

0.53 

65 

279.9  5.23 

1.87 

75 

270.1  6.47 

2.40 

80 

264.2  7.14 

2.70 

85 

256,9  4.29 

1.67 

» 
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We  simply  write  down  the  analogous  results  for  T2(0q): 


V2^®o^  =  dSj(cosaj)rj(l-r2) 


where 


=  S^ta2(e„)  t 


°  l<ISj(cosaj)  rj{l-r2) 


(3.11) 


(3.12) 


and  the  residue  term,  A2(eQ),  contains  higher-ordsr  absorption  and  scatter 
processes.  The  subset  has  been  further  subdivided, 
such  that  the  ray  k2j  either  clears  the  surface  [(Jlj)^]  or  strikes  the  surface 
again  I(Qp^].  The  quantity  T2(0q)  is  the  emission  plus  reflection  term^  in 
the  lower-bound  calculation  of  rough-surface  emissivity  and,  again,  a  computer 
program  is  available  for  its  evaluation  (Tables  I  and  II).  The  quantity 
a2(0Q)  is  a  contribution  to  the  absorptivity,  by  construction,  and  it  can  be 
calculated  numerically.  Thus  A2(0q)  can  be  calculated  from  Eq.  (3.11).  We 
neglect  the  remaining  terms  T^,  ...  T^j^  relative  to  Tj  and  T2;  this  is  a  valid 
approximation  for  the  rms  surface  si  opes .appropriate  to  the  ocean  surface. 

In  summary,  we  note  that  Tj,T2,  •••  can  be  thought  of  as  channels 
through  which  the  incident  energy  is  funneled  into  an  infinite  number  of  pro¬ 
cesses.  For  example,  the  double  transmission  process  of  Fig.  4  is  funded  by 
the  energy  budget  associated  with  Tj.  It  is  not  necessary  to  evaluate  the 
magnitude  of  these  higher-order  processes  on  an  individual  basis,  but,  instead, 
it  is  meaningful  to  lump  them  together  in  the  A..  For  example,  the  fraction 
of  the  incident  energy  funneled  into  channels  Tj  and  T2  can  be  calculated,  so 
Aj  and  can  be  determined.  We  will  simply  assume  here  that  the  higher-order 
absorption  processes  in  Aj,  A2  (and  T^,  ...  T^^^)  are  negligible  compared  to  the 
smallest  of  (ai>2iR>a2^  slopes  associated  with  the  ocean  surface. 

Then,  a  meaningful  expansion  for  the  absorptivity  is 

3(0^)  =  ai(9o^  ^‘1R‘®0^  ®2^®o^ 
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with  (aj,ajp,a2)  defined  formally  by  Eqs.  (3.8),  (3.9),  and  (3.12)  respectively. 

In  passing,  we  note  that  every  term  in  the  T  equation  is  positive.  The 
technique  of  this  section  was  to  further  subdivide  each  of  the  into  a  series 
of  positive  terms,  some  of  the  elements  to  be  grouped  under  a(0Q)  and  the  rest 
under  6(0^)  [Eq. (3.4)1.  The  net  result  is  that  each  term  on  the  right-hand  side 
of  Eq.  (3.13)  is  positive.  Thus,  by  taking  just  the  first  m  terms  as  an  approxi¬ 
mation,  we  have  a  lower  bound  to  the  absorptivity. 

IV.  PARTIAL-SURFACE  REPRESENTATION  FOR  THE  TRANSMISSION  COEFFICIENT 

The  absorptivity  provides  no  knowledge  of  the  angular  distribution  of 
transmitted  energy.  We  define  the  transmissivity,dJl  t(k,kQ;£,)  as  the  fraction 
of  energy  incident  from  ]<  which  is  refracted  into  solid  angle  dfl  about  the  ob- 
server  direction  k,  for  an  observer  at  depth  £  (measured  from  the  mean  surface 
height).  It  is  conventional  to  deal  simply  with  t,  the  transmission  coefficient. 

As  with  the  absorptivity,  there  is  a  partial -surface  expansion  for  t.  Indeed, 
as  a  consequence  of  ray  theory,  the  partial  transmission  coefficients,  tj,  t^j^, 
t2  are,  at  least  formally,  simple  extensions  of  the  corresponding  a^,  aj|^,  a2 
of  Section  III. 

We  will  first  detail  the  derivation  of  the  direct  transmission  coefficient, 
tj,  for  a  one  dimensional  rough  surface.  The  appropriate  partial  surface  is 
ny*,  i.e.,  every  surface  element  dSj  in  Qj  is  visible  to  the  source  and  to  the 
observer.  The  pow-r  incident  on  dSj  is  (cosaj)dSj,  and  the  fraction  transmitted 
is  dSj(cosaj)  [l-rj(cosaj) ].  Since  dSjeJlj,  every  refracted  ray  clears  the  (under) 
surface.  We  account  for  the  absorption  along  the  propagation  path  of  the  ray  by 
a  multiplicative  exponential  absorption  factor,  exp[-K(£+Zj)sec0jj] ,  where  K  is 
the  absorption  coefficient,  Zj  is  the  surface  height  of  point  1,  and  0jj  is  the 
angle  of  propagation  for  the  refracted  ray.  We  must  also  multiply  by  the  proba¬ 
bility  that  the  ray  reaches  the  observer  stationed  at  angle  0;  because  of  the 
ray  optics  and  the  assumed  homogeneity  of  the  medium,  the  probability  density 
is  simply  a  delta  function  6(0-0^^}.  Finally,  we  sum  over  all  dS^efly: 

Soti(0,0o)  =  |dSj(cosaj)(l-ry(cosaj)]exp[-K()l+Zj)sec0yy]6(0-0jy)  ,  (4.1) 

fiy 

where  S^,  defined  in  Eq.(3.3),  is  included  on  the  left-hand  side  for  normalization. 
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'  ty"'  -^r^ ,  --Ji^fa^^ 


l| 

I 

i 

► 

I 

I: 
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The  analogous  expression  for  the  direct  transmission  coefficient  applica&ie 
to  a  general,  two  dimensional,  random  rough  surface  is 

Sj^tj'^(k,ko)  =  dSj(cosaj)(T^'^)expl-K(Jl+2j)secejjl6(0“0iT)5('l>-<l>ij)  •  (4.2) 

Jlj 

Here,  the  superscripts  reflect  the  polarization  character  of  the  problem.  Speci¬ 
fically,  the  quantity  t^'^(l<,!<  )  determines  the  fraction  of  the  incident  radiation 
from  Kq,  in  polarization  state  i,  which  is  transmitted  into  an  infinitesimal  solid 
angle  about  observation  direction  t,  with  final  polarization  state  j.  The  fraction 
of  incident  power,  with  polarization  state  i,  transmitted  across  an  arbitrarily 
oriented  plane  surface  into  final  polarization  state  j  is  represented  by  for 

•  • 

cylindrical  symmetry  the  incident  polarization  is  retained  after  refraction  and 
reduces  to  (l-r^)fi^j.  The  general  form  of  is  derived  in  Appendix  B.  Equation 
(4.2)  will  usually  be  the  major  component  of  the  exact  transmission  coefficient. ' 

We  will  simply  state  the  results  for  t^j^  and  t2,  and  we  will  restrict 
consideration  to  one  dimensional  roughness: 

» 

SgtiplO.Oo)  =  dSj(cosaj)(l-rj)exp[-K{Zj-Z2)sec0j.p] 


‘TR 


x  r2  exp{-K(£+Z2)sec02p] 


X  5(0-021^)  , 


(4.3) 


^o^2^®’®o^  =  dSj(cosaj)rj(l-r2)  exp[-K(£.+Z2)sec02j] 

iL> , 

X  6(0-02j)  (4.4) 


(«;)t 


The  series  for  the  transmission  coefficient  has  a  form  analogous  to  Eq,  (3.13): 

t(0,eo)  =  tj(9,0g)  +  tjj^(9,9Q)  +  t2(0,9o)  +  ...  (4.5) 

Each  term  in  Eq.  (4.5)  is  positive,  so  the  retention  of  the  first  m  terms  pro¬ 
vides  a  lower  bound  to  t{QyQ^). 
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V.  DIRECT  TRANSMISSION  COEFFICIENT  FOR  RANDOM  ROUGH  SURFACES 

The  formalism  of  Section  IV  i'  especially  suitable  for  application  to  a 
random  rough  surface.  Thus,  we  will  view  c{x,y)  as  one  of  an  ensemble  of  possible 
surfaces  generated  by  a  stationary  random  process  and  calculate  the  ensemble 
average  of  the  transmission  coefficient.  The  averaging  process  replaces  the  deter¬ 
ministic  surface  classes  by  known  probability  functions,  and  it  also  provides  for 
an  integration  over  the  delta- function  factors. 

We  will  consider  only  the  direct  transmission  contribution  t^^ 
this  section.  We  can  convert  the  integration  over  Qj  into  an  integration  over 
the  entire  surface  I  by  defining  two  functions  which  take  on  only  the  values  of 
zero  and  unity.  Thus,  we  define  an  illumination  function  e(x,y;k^),  which  has 
value  unity  If  the  surface  point  at  (x.y)  is  illuminated  from  direction  and 
which  has  value  zero  if  that  point  is  shadowed.  In  addition,  we  define 
v(x,y;kjj)  to  have  value  unity  if  the  transmitted  ray  which  leaves  point  (x,y)  in 
direction  k^j  intersects  the  surface  at  some  other  point,  whereas  it  has  value 
zero  if  k^j  does  not  intersect  the  surface  again.  With  these  definitions, 

Eq.  (4.2)  becomes 


Sot|'’(k,ico)  =  |dSj(cosaj)T^^  exp[-K()l+Zj)seceij3 
Z 

X  e(x,y;kjj)ll-v(x,y;kjj)] 

X  6(e-0j-j-)  6(())-it>jy)  .  (5.1) 

Now,  the  coordinates  of  the  transmitted  ray  for  surface  element  dSj  are 

^®1T’'^1T^’  Snell's  law  we  have  e^y  =  0iT(®o’*lx»^ly^ 

(|)1T  =  <J>iT^®o’^lx’^lv^*  ^hx’^ly^  components  of  surface  slope  at 

point  1.  The  surface  element  is  dSj  =  (l+s^j^  +  Sjy)*^  dx^dy^,  and  cosa^  and 
depend  only  on  (Xj,yj)  through  the  slops  (Si^’S^y).  Therefore,  the  integrand  in 
Eq.  (5.1)  depends  on  (xj.y^)  only  through  Zy(xy,yj), 

e(Xpyj),  v(xpyj).  For  any  stationary  random  process,  there  are  no  preferred 
points,  50  the  appropriate  probability  distribution  is  independent  of  (Xj.y^. 
Thus,  the  averaging  procedure  leaves  a  trivial  coordinate  integration: 


So<‘P(^4>)  ■  I  I 

•[_  .CO  ^co  -00  IsQ  X 

X  ei(l-Vi)fi'*(Si,.Siy) 

X  expt-K()l+Zj)secejy)] 

X  6(0-6j^)6(<>-(|.^y) 

»  (2L)(2L')  jdZjjdSjj^pSjy  PtZi»Six»Siy.ei(ko)  »  l.Vj{ICjj)  =  Ol 

MQC  -00  >-00 

X  exp[-K{fc+Zj^)sec0jj]6(e-0jj)6((^-^jj)  , 

(5.2) 


where 

'l'’('lx'*ly>  '  (l'"®lx'"*ly>''  <““l'  ^''’<*lx>®ly>  ’ 

The  segments  along  the  x  and  y  axes  cut  out  by  the  incident  beam  are  taken 
here  as  2L  and  2L',  respectively,  so  that,  in  these  terms,  =  (2L)(2L')cos0q). 
Also,  we  note  that  the  ray  kjj  must  intersect  the  surface  if  0jj  >  -n/Z  (Fig.  2). 
For  <  Tr/2,  the  probability  that  the  ray  kjj  does  not  intersect  the  surface 
elsewhere  is  equivalent  mathematically  to  the  probability  that  the  point  (x^y^) 
is  illuminated  (by  an  imaginary  source  below  the  surface)  from  the  direction 

On  the  basis  of  these  remarks,  Eq.  (5,2)  can  be  expressed  in  the  following 

form: 
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Mh 

t/f 


(t\hU,)) 


(sece„)pzjpsj^|dsjj,  PtzpSj^.Sjj,.  ei(k(,.-k„) 


=  11 


-CO  «00  -00 


X  U(ti/2-6jj) 

X  exp[-K(Jl+Zj^)secej-|.] 
X  6(6“62j)6{<)>“<j)2j)  i 


(5.4) 


where  u  is  a  step  function  with  value  unity  for  positive  argument  and  value  zero 
for  negative  argument. 

The  above  multiple  integral  is  quite  complicated  because  of  the  integration 
over  height  z^.  It  is  sensible  to  approximate  the  factor  Zj  in  the  exponential 
argument  by  the  average  value  of  the  "illuminated"  surface  height,  S,j,  defined  by 


(5.5) 


where  the  probability  density  of  heights  is  conditional  on  the  illumination  of 
point  1  from  both  the  and  -kjy  directions.  The  gross  effect  of  absorption  is 
therefore  retained  while  the  analysis  is  greatly  simplified.  The  presence  of 
the  delta  function  permits  us  to  replace  e^j  by  0,  so 


^tj^(t,k^)^  =  (sec0Q)exp[-K{it+2’j|^)sac0] 

*  P"lf=lxPhy'’'^l-Slx-hy-n<''o'-'‘ 


IT 


)  =  1] 


—00  -CO  —00 


X  fl^(six’Siy)«(TT/2-ej.|.) 


X  <S(0“6j^j)6(<j>”<l*jj) 
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(tj^(k.ko)>  =  {sec0Q)exp[-K(£+x-j^)sece} 

••00  >>00 

x  fp{Sj^.Sjy)6{0  ejT)6{(}>-<|>iy)  ,  (5.6) 

where  the  last  step  follows  from  the  assumed  normalization  of  the  probability 
density. 

The  remaining  integrations  can  be  carried  out  by  a  change  of  variables  from 


(Slx.Sly)  to  (0j-j-,<|)lj). 

The  equations  of  transformation  are. 

from  Appendix  C, 

(sin0Q-  n  sin0j.|.cos(i)j^j) 

hx 

(n  cos0^j-  cos0q)  ’ 

-  n  sinBjj  sin^jj 

(5.7) 

^ly 

(n  COS0JJ-  cos9q)  * 

and  the  Jacobian  of  the 

transformation  is 

J  (Ojjj't’jy)  ~ 

O  A  /N 

n  sin0jj[n-{kjj* 

(5.8) 

o  • 

(n  cos0j.j.-cos0q) 

Here,  n  is  the  real  part  of  the  index  of  refraction.  One  must  be  very  careful 
at  this  stage  because  the  integration  limits  on  the  (Qjtj'I'ij)  integrals  are  not 
simp!'/  0  <  9i-r  <  tt/2,  0  <  (j)jj  <  Zv.  This  follows  because  the  orientation  of  a 
plane  surface  which  connects  the  desired  direction  for  transmission  with  the 
given  direction  of  incidence  may  not  be  allow'^d  on  physical  grounds.  More 
specifically,  '"are  are  two  conditions  on  which^follow  from  physical 

constraints  on  this  problem.  First,  we  require  that  n^*  z  >  0  everywhere,  i.e., 
the  angle  between  any  surface  normal  and  the  z  direction  must  be  less  than  90°. 
This  restriction  arises  because  we  are  assuming  a  single-valued  normally  dis¬ 
tributed  surface,  a  property  of  which  is  finite  slopes  everywhere.  But,  from 
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)  Appendix  C,  we  see  that  this  condition  implies 

cose. 

cos0j^  >  — —  ,  (5.9) 

A  second  condition  is  the  illumination  requirement  (-k^*  n^)  >  0;  i.e.,  a  surface 
point  cannot  be  illuminated  by  the  incident  ray  unless  the  angle  between  the 
normal  to  the  point  and  the  direction  of  illumination  is  less  than  90°.  This 
obvious  condition  imposes  a  very  severe  restriction  on  (6ji-,!|)jj) : 

A  1 

''o**'lT^n  .  (5.10) 

The  simple  way  to  incorporate  these  two  restrictions-on  the  transmitted  direction 

/  A  A  2  ^ 

is  by  use  of  two  more  step  functions,  ulcose^^ - — |  and  u(kj  kiT--);  i.e., 

fCo  M  ,t!/2  ,  cos9  . 

j^^lx  ^^ly  ^  *^^®1T*'^1T^  ujcosGjj  | 

-00  -00  0  0 

■is  X  u(k„-  ■  n|  •  (5.11) 

The  step  functions  automatically  impose  the  proper  integration  ranges. 

The  application  of  Eqs.(5.7),  (5.8),  and  (5.11)  to  Eq.(5.6)  yields 

t|^(k,kjj)  =  (secSp)  exp(-K(Jl+i.2)sec0] 

'Tr/2  *2Tr  a  a  ^ 

^  ^^IT  ^*^1T  u(ncos0jj”  cos0Q)u(kp*  ”) 

J  J 

0  0 

^  fj'^  5{9“92j)5(<l>”<j>2-j-) » 
which  simplifies  immediately  to 

I 
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(ti^(k,ko)),  -  (sece^)  {expt-K(ji+]lj)sece]}fp(s°j^,s°y) 


X  u(ncos0-cose„)u(k  •  k- 4-) 
0  0  n 


n^[n-(k  •  k)] 


(ncosS-cose^)' 


■3 


(5.12) 


where  we  have  used  the  properties  of  the  delta  function  and  the  normalization 
r7r/2  f2ir 


fir/c  rCTT 

Jd0iTjd4.1T 


sine^T  6(0“0iT)5(<|f~<t>iT)  =  1 


IT 


IT' 


Geometrical  quantities  in  Eq.(5.12)  which  need  to  be  identified  explicitly  are 

„o 


(sin0Q-nsin0cos4)) 


’lx  "  (ncosG-cosOg) 

.0  _  ■nsin0sln4>  • 
*ly  "  (ncos0-cose^) 


'  and 


-  ( 

whence,  after  some  algebra, 

‘lx  ^ly^ 


.2. 


1^ 

,  .TV  ^ 


o2_  qZ^H  (n  +,l-2n(kg«.k)] 
(1+Si„+Si„)  -  (ncos0-cos0Q)  “ 


cosa 


ln^+l-2r,(kQ-k)]^" 


(5.13), 


(5.1,4) 


(5.15) 

(5.16) 


oriented 
and  cosa 


From  Appendix  B,1  the  fraction  of  energy  trifnsmitted  acros's  the  properly 

c®  ^  -ic  a  nn  nf  a®  h®  +0 


plane  .surface,  ^  '^(Sj^.s^y),  is  a  function  of  a^,  bp  Cp  t^j,  ty,  cosaj 
l(Sp,Siy).  The  following  list  specifies  each  of  these  quantities  in 


terms  of  known  factors: 
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I 


I/O  0  \ 

COSaj(s^j^,Sjy)  = 


t{n  -l)+cos'^ajl 


(5.17) 


r\  r\ 

Ik.xkl 


A,  A 

-k-v„ 


;  ii„  =  y  . 


1^1  :  v„  = -cose^  X  +  sine^z  , 


-(I<,X  k)‘{i<  X  i) 
(lic^jX  kl)(|k  X  z|) 

-(k^x  k)* [z-(ic»z)k] 
(ll<  X  (c|)(|k  X  z|) 


(5.18) 


(5.19) 


(5.20) 


(5.21) 


tj  =  1  +  R^(cosa5) 


^0  _  l+Ry(cosaj) 

ty  =  ^ 


The  probability  density, has  the  form 


^^^Ix’^ly*  =  11  =  p(s5j^,s°y)P[ej(ltp,-t) 


(5.22) 


(5.23) 


IjsJ^.sJy]  .  (5.24) 


where,  for  the  principal  directions,  the  probability  density  of  slopes  is 

f 


(5.25) 


and,  from  Eq.  (6.11)  of  the  following  section,  the  probability  of  illumination 
conditional  on  the  slope  value  is 


P[ei(ko,-i<)  -  ll(S]x*Siy)l  -  fi+2(B  +b)3  r[i+2(B +B)] 


r(l+2Bjj)r(l+2B) 


(5.26) 


/ 


The  quantities  r(x)  are  Gamma  functions  of  argument  x,  and 


B  =  [4(tt)^*V]  \exp(-V^)  -  V(l-erf  V)]  ,  (5.27) 

with 

V  =  - ^ - j:  (5.28) 

C2((c^cosV(c^sin^<t.)l  ^ 


The  quantity  follows  from  B  by  replacement  of  (0,4))  by  (0o»‘J’o  “  Equations 
(5.13)  -  (5.28)  specify  all  quantities  in  Eq.  (5.12)  except  the  average  il¬ 
luminated  height  in  the  exponential  absorption  factor.  This  factor  is  easily 
derived  from  the  shadowing  theory  of  Section  VI: 


S'l(k.koy  =  (77/2)  ^0  ^+8^+5) 


(5.29) 


where  o  is  the  rms  height  for  the  surface. 

The  numerous  factors  which  we  have  just  listed  must  not  be  allowed  to 
obscure  the  overall  simplicity  or  Eq.  (5.12).  Thus,  there  are  no  integrations 
to  carry  out;  instead,  the  direct  transmission  is  simply  proportional  to  the 
probability  density  of  that  particular  slope  which  connects  k  and  k^  by  Snell's 
lav  The  probability  that  this  characteristic  slope  is  visible  to  both  source 
and  observer  also  enters.  An  additional  feature  of  interest  is  the  product  of 
the  two  step  functions,  which  can  yield  zero  for  some  combinations  of  t  and 
Thus,  it  is  not  always  possible  to  connect  a  direction  of  transmission  with  a 
given  direction  of  incidence  by  any  physically  allowed  plane-surface  element. 

This  is  most  easily  demonstrated  for  the  special  case  of  the  one  dimensional 
rough  surface.  Then,  the  allowed  range  for  the  transmission  angle  8  is  easily 
shown  to  be 

0^  -  [7T/2-sin"^(^)]  <  0  <  7T/2-sin"^|^^^j  .  (5.30) 

For  example,  if  the  angle  of  incidence  is  0^  =  20°,  then  there  is  no  direct 
transmission  for  0  <21°  and  0  >  45°.  Of  course,  the  higher-order  processes, 
tji^  and  t2,  will  contribute  something  to  these  "null"  ranges. 
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VI.  SHADOWING  PROBABILITIES 


An  element  of  surface  has  only  a  certain  probability  of  being  illuminated 
by  radiation  incident  from  directions  other  than  normal.  For  example,  the  angle 
between  the  surface  normal  and  the  direction  of  incidence  may  exceed  90°,  so 
that  the  slope  of  the  element  totally  rules  out  illumination.  Or,  the  topography 
of  the  intervening  surface  may  simply  cast  its  shadow  on  the  test  element.  The 
probability  that  a  point  on  a  surface  is  illuminated  depends  on  the  angle  of  in¬ 
cidence  as  well  as  the  height  and  slope  of  the  surface  point. 


In  its  application  to  scattering  theory,  the  appropriate  shadowing  formal- 
ism  was  found  to  depend  on  both  the  angles  of  incidence  and  observation.  The 
application  to  transmission  theory  will  be  bistatic  as  well.  The  shadowing 
theory  for  transmission  differs  substantially  from  that  for  scattering,  however. 
Thus,  a  surface  high  point  is  more  likely  to  be  illuminated  by  the  incident  ray 
than  a  surface  low  point,  but  the  surface  high  point  is  less  likely  to  be  visible 
to  an  arbitrarily  oriented  observer  below  the  surface  than  the  corresponding  low 
point  (Fig.  1). 


We  consider  the  incident  ray  located  in  the  x-z  plane,  and  the  arbitrary 

->.0  y 

observer  direction  k^j,  defined  as  in  Fig.  2.  From  Sancer,  the  probability  that 
surface  point  I  is  illuminated  from  k^,  given  its  height  and  slope,  is 


S(9o|Zi.Six-Sly) 


1+erf 


(6.1) 


where 


Bq  =  (4(Tr)\)'^[exp(-v2)  V^)] 


COt0Q 


(6.2) 

(6.3) 


Again,  the  quantity  u  in  Eq.  (6.1)  is  a  step  function,  which  is  unity  il  the 
argument  is  positive  and  which  is  zero  if  the  argument  is  negative.  The  qae'tity 
B^  varies  from  zero  to  infinity  as  6^  varies  between  normal  and  grazing  lAeidertce 
Also,  the  height  dependent  factor  in  Eq.  (6.1)  varies  from  zero  to  unity  (for 
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nonnormal  incidence)  as  varies  from  minus  to  plus  infinity.  This  is  the 
desired  behavior  for  the  incident  direction,  but  for  the  transmitted  (observed) 
direction  we  want  the  inverse  variation;  i.e.,  as  varies  from  plus  to  minus 
infinity,  the  probability  of  illumination  should  increase  from  zero  to  unity. 

By  use  of  the  same  analysis  applied  to  the  derivation  of  Eq.  (6.1),  we  conclude 
that 


S(0jt»'I>itIzi.Sjx.Si  )  =  u(-kj.p*n^) 


1-erf 


■1 


p  h 

/j 


2B 


IT 


(6.4) 


0 


and  Bjj  follows  from  B^  by  replacement  o"  by  where 

cote. 


'IT 


'IT 


2((Cx)co^^<<>1T 


(6.5) 


The  desired  quantity  is  the  bi static  probability  of  illumination.  Fortu- 
nately,  because  and  k^j  are  associated  with  different  quadrants  of  space,  we 
can  assume  statistical  independence: 


Cl 


5(^0, kiflzi, Six, Sly)  =  5^®ol^l’^lx’^ly^^^®lT’‘^’lTlh’^lx*^ly^ 


>\  /s  /s 

=  u(-RQ*ni)u(-kiT*ni) 


1+erf 


p  ^  I 
(20^) 


2B.r 


1-erf 


■1 


p  h 

(2o^) 


2B 


IT 


(6.6) 


Equation  (6.6)  provides  more  information  than  we  need,  as  the  transmission  co¬ 
efficient  formalism  of  Section  V  requires  only  the  probability  of  illunination, 
conditional  bn  the  slope  value.  Therefore,  we  simply  average  Eq.  (6.6)  over 
heights.  For  a  normally  distributed  surface,  the  expression  is 


S(ko.kjTisix,Siy)  =  JdhP^h^S(k^,ki^lzj,Sjx,Siy) 


(6.7) 
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By  means  of  the  substitution 


1 

y  =  7 


and  the  property 


d/dz  erf 


7  ^ 

l-(2<i^)  . 


\ 

1+erf  .  ’ 

^  (2a2)“ 


2  P(z)  . 


we  have 


1  2B.  2B 


AA  AaT^ 

S(ko,kiTls^x*Sly)  =  u(-kQ*n^)u(-kj^‘nj)Jdy  y  °(l-y) 


IT 


(6.8) 


u(-kQ*nj)u(-kj^*nj)  r{l-f2B|j)r(l-i-2Bj.|-) 
[l+2(Bp+Bj.f)]  r[l+2(Bjj+Bj^)] 


(6.9) 


The  quantities  r  are  Gamma  functions/  and  the  integral  in  Eq.  (6.8)  is,  of 
course,  the  standard  representation  for  the  Beta  function  (BQ,Bjy  ^  0). 

It  needs  to  be  emphasized  here  that  Eq.  (6.9)  is  valid  for  any  arbitrary 
kjj.  In  the  application  to  the  transmission  problem,  however,  k^y  is  not  an 
independent  vector,  but  it  is  related  to  k^  and  tne  surface  normal  n^  by 
Eq.  (C.4);  i.e. , 


n(-kjj‘nj)  =  (-kp*nj)  +  {[n  -1+cos  a^l  -cosoj} 


(6.10) 


It  is  easily  seen  from  this  relation  that  if  (-k^^nj)  >  0,  then  (-kjj*nj)  >  0  as 
well.  Therefore,  for  the  transmission  problem,  we  can  set  ‘^("kjj'n^)  equal  to  unity 
without  loss  of  generality.  In  addition,  we  have  already  incorporated  the  step  func¬ 
tion  u(-kQ'n^)  in  the  transmission  theory  of  Section  V  in  the  context  of  the  change 
of  variables  from  (s^j^.s^^)  to  (92j,<})2j).  Thus,  the  slope  coordinates  (s2x»S2y) 
in  Eq.  (5.12)  can  be  assumed  to  be  physically  allowed.  Therefore,  we  now  convert 
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Eq.  (6.9)  to  the  notation  of  Section  V  in  order  to  obtain 


PU 


£1(^0, -k)  = 


0  „o 


r{l+2BQ}r(l+2B) 


11+2(Bq+B)]  rIl+2(Bjj+B)] 


(6.11) 


V!I.  REFLECTION  FOLLOWED  BY  TRANSMISSION 

Radiation  incident  on  a  rough  surface  will  undergo  multiple  scatter  due  to 
that  roughness.  At  each  point  of  contact  between  ray  and  surface,  there  will  be 
a  contribution  to  the  transmission  of  energy  across  the  boundary.  The  coefficient 
yt,)  of  Section  V  accounts  only  for  the  transmission  into  the  medium  at  the 
initial  point  of  contact  [the  ray  in  Fig.  (3)].  In  this  section,  we  will 
write  down  the  coefficient  ^t2)  appropriate  to  transmission  at  the  secono  point 
of  contact  [the  ray  k2j  in  Fig.  (3)1.  Only  multiple-scatter  ray  diagrams  of 
order  two  and  higher  will  contribute  to  this  term.  A  good  estimate  of  the  im 
portance  of  this  process  can  be  obtained  from  Section  III  where  a  numerical 
evaluation  of  the  quantities  Tj(0q)  and  T2(0q)  has  been  carried  through.  From 
Tables  I  and  II,  we  see  that  the  ratio  (T2/Tj)  has  a  broad  maximum  of  roughly 
3%  at  the  large  angles  of  incidence. 

The  partial -surface  representation  of  t2(0,0p)  is  given  by  Eq.  (4.4). 

Because  of  the  cumbersome  mathematics  involved,  we  will  restrict  the  details  to 
the  one  dimensional  rough  surface,  c(x,y)  =  c(x).  Also,  the  averaging  process 
is  greatly  simplified  if  we  assume  points  1  and  2  to  be  statistically  inde¬ 
pendent.  This  appears  to  be  an  excellent  approximation  because  the  slopes  at 
points  1  and  2  will  usually  be  appreciably  different  in  value,  the  characteristic 
of  a  separation  distance  larger  than  the  correlation  length.  By  use  of  techniques 
similar  to  those  employed  in  Section  V,  we  have 
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(t2(0,9Q),)  =  (seceQ)|dSjpS2  PlSpS2,Ii(eQ)  ,1^2(61) 

•00  *00 

X  u{Ti/2-|e2j|  )Ti'2^^1’^2^ 

X  exp  {  “K[A+Jl2(6i)^sec0  } 
X  6  [6“621"(  0q»s  j  >52)  1  > 


where 


(7,1) 


0  • 

f2(Sj,S2)  =  (1+sp  (cosaj)rj^lcosaj(sp]{l-r2fcosa2(Sj.S2)l}  ,  (7.2) 


£2(61)  =  pZ2  7.2ptz2|li{0o),  112(81). * 


(7.3) 


The  quantity  P  is  a  joint  probability  density  on  slope  and  illumination  for  each 
of  the  two  scatter  points.  The  symbol  112(61)  denotes  the  likelihood  that  the 
intermediate  ray  ki  actually  intersects  the  surface  again  at  point  2;  it  must  be 
included  as  a  fifth  random  variable  in  any  consistent  theory  for  double  scatter. 

As  the  argument  of  the  delta  function  depends  on  02j  rather  than  S2.  it  is 
clearly  advisable  to  change  variables  from  (Si,S2)  to  (Si,e2j).  The  equation  of 
transformation  can  be  deduced  from  Appendix  C: 


(nsin02i-+sin0j^) 

(ncos02j+cos0i) 


(7.4) 


where 


9,  =  -0„-2tdn  s,  , 

1  0  1 


(7.5) 


or,  inversely, 


|6n+6,  , 

-tanj-^w — ) 


(7.6) 
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The  Jacobian  of  the  transformation  is 


J(si) 


n[n+cos(02j'’6j)l 

(cos6j+ncos02j)^ 


(7.7) 


a  function  of  02j  and,  implicitly,  of  s^.  As  was  the  case  with  the  direct  trans¬ 
mission  coefficient  (Section  V),  the  -ange  of  integration  for  02j  is  restricted 
by  the  illumination  condition,  (-kj*n2)  >  0,  and  by  the  condition  (n2*z)  >  0, 
the  latter  a  consequence  of  single-valued  normally  distributed  surfaces.  An 
additional  complication  arises  because  some  values  of  Sj  lead  only  to  values  of 
|02tI  larger  than  Tr/2  (e.g.,  the  vector  k^j  in  Fig.  3).  The  following  expression 
for  ^t2)  accounts  for  the  change  of  variables,  the  resulting  integration  over  02y, 
and  the  imposition  of  the  preceding  restrictions: 


(t2(e,6„)> 


(sec0Q 


jdSj  J(s^)u(-Xj-6)F(Sj,S2) 
sin"^(i)  <  0j  <  Tr/2 


dSj  O(Sj)u(0+Xj)u(-Xj-0)F(SpS2) 


0j>Tr/2 


+  jdSj  J(Sj)u(e-x^)F(Sj,S2) 
-Tr/2<0j<-sin”^(^) 

+  ds,  O(Sj)u(0-Xj)u(xi-0)F(SpS2) 
0j<-Tr/2 


(7.8) 


where  the  integration  ranges  for  the  slope  integrals  are  given  in  terms  of  equi 
valent  ranges  for  the  angle  0^  (Eq.  7.6).  The  quantities  x|[.  in  the  step- 
function  arguments  are  defined  as 
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u  -i/coseii 

xj  =  7r/2  +  sin 

(7.9) 

Xi  =  {Tr-|9i|)  -  [ir/2-sin  ^i)] 

(7. 1C) 

The  quantity  FCs^.Sg)  is 

F(sj,s°)  =  PtSj,s°,Ij(0Q).  Ii2(6i).l2(®^^^2^^1’^2^ 

X  exp{-K[£,+Jl2(0j)^sec0} 

(7.11) 

vhere 

p  (sin0j+ns1n0) 

5  —  —  '  » 

(cos0^+ncos0) 

(7.12) 

and  f2  is  given  by  Eq.  (7.2),  evaluated  at  $2  =  Sg.  The  arguments  of  the 
flection  coefficients  can  be  expressed  most  usefully  as 

re- 

2 

cosaj(sj)  =  (1+sp  (cose^j  -SjSin0jj) 

(7.13) 

cosa2(Sj,  ,S2)  =  -[l+(sp^]  (cosO^-SgSine^) 

(7.14) 

The  fifth-order  probability  density  can  be  written  as  the  product  of  five 
probabi 1 i ty  f uncti ons : 

P[s j jS2»I2^(cq) j  I ~  *  (7*15) 

where 


1  =  Ptlj(0p)|Sj]  , 

(7.16) 

(7.17) 
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P3  =  p[s°|sj,ii{eQ).  112(61)^ 


(7.18) 


P4  =  p[i2(e)lspS°.ij{e^).ij2(0i)^  • 


(7.19) 


The  probability  density  P(s2)  is  the  probability  density  of  slopes  for  a  normally 
distributed  surface.  The  mathematical  relations  for  PpP2»  and  P3  can  be  found 
in  the  literature.^  The  remaining  P^  asks  for  the  probability  that  an  observer 
in  the  medium,  oriented  at  ang'e  &,  receives  the  refracted  ray,  given  the  values 
of  slope  at  points  1  and  2  and  given  the  facts  of  illumination  specified  by 
ll(0o)  and  Ij2(0i)*  The  analytical  results  for  P^  are 


(1+2B) 


Qji  <  tt/2 


(l+2B 


j  r(i+2Bj)r(i+2B) 
il-:-2(Bj+B)]  r[l+2  (Bj+B)] 


ISjl  >  Tr/2  ,  (7.20) 


where  the  B's  have  the  functional  dependence  given  by  Eq.  (5.27)  and  V,Vt  are 
given  by 


% 


(7.21) 


I cot(n-62) I 


(7.22) 


The  only  remaining  quantity  to  be  specified  is  il2(ei).  the  average  illum¬ 
inated  height  of  point  2.  The  results  are 
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^  ^  mi^l.llJjM^,JUP.i!!!l  jyiJ,kUi^iiW  w.;iiBi!.w^i  >  ^ 


•5  ' 


il2{9i) 


_J±_  (Tr/2)^*  0 
(1+B) 


le^l  <  V2 


,  |Tld‘''2«Bi-2BH  r(H2Bj)r{2+2B) 

=  - - - i - = 

[r(2'5-2Bj+2B)  -  r(l+2Bj)r{2+2B) 


Bjl  >  tt/2  .  (7,23) 


VIIL  TRANSMISSION  FOLLOWLD  BY  REFLECTION 


The  partial  transmission  coefficients  tj  and  t2  describe  a  transmitted  ray 
which  travels  directly  to  a  hypothetical  observer;  i.e.,  the  transmitted  ray  does 
not  intersect  the  surface  again.  However,  Fig.  4  demonstrates  a  higher-order 
process  with  transmission  at  point  1  followed  by  a  surface  reflection  at  point  2 
(the  ray  i<2R^’  section  we  derive  the  average  transmission  coefficient 

associated  with  this  process.  This  effect  is  included  only  for  completeness 
as  its  contribution  to  the  transmitted  radiation  pattern  is  expected  to  be  quite 
small.  This  follows  because,  as  Fig.  4  implies,  the  angle  9jj(0jj  <  Tr/2)  must  be 
quite  large  f  the  ray  is  to  have  much  chance  of  intersecting  the  surface  again. 

But  it  is  easy  to  show  that  a  large  value  of  requires  a  combination  of  large 
angle  of  incidence  and  an  element  of  surface  with  steep  slope.  Thus,  (^r)  can 
contribute  significantly  only  for  near-grazing  incidence  and  the  very  roughest 
of  ocean  surfaces.  We  restrict  the  details  in  the  following  to  the  one  dimensional 
rougn  surface,  ?(x,y)  =  ^(x). 


The  partial -surface  representation  for  tj|^(0,0Q)  is  given  by  Eq.  (4.3). 

The  assumption  that  points  1  and  2  are  statistically  indepet  dent  again  appears  to 
be  an  excellent  one,  for  the  slope  values  at  the  two  points  must  be  very  different. 
This  implies  a  separation  distance  gi  ater  than  the  correlation  length.  The 
average  value  of  t,(^  follows  by  the  techniques  of  Section  V: 


f 


~  ds^  ds^  P Is «S2* I j(6q)  ,I ^^2(0 j^j)  ^  ^ 


X  u(Tr/2- 1 02R|)fiR(Si  ,82) 

'  \ 

1 

1  1 

X  exp{-K[A3(0ii.)-.<l4(0ij)  ]sed0iy} 

'  ' 

' 

X  expt-K(Jl+il4)sec0] 

/  ’  ‘ 

(8.1) 

where 

1 

flR(Si,S2) 

2  ^ 

=  (l+sp  (costti)  [l-ri(cosoti)3r2(kij*n2)  , 

(8.2) 

I  ,  ' 

•  ^  1 

=  .  pZi  Zi  Ptzi|Ii(0Q).,  112(^17)^  > 

^mCO 

'  (8.3) 

2.4(611) 

I 

#UO  1 

=,  jdZg  Z2  Ptz^l  112(617).  * 

-09 

(8.4) 

and 

'  •  .  .  *  ' 

i 

Cosai, 

=  (l+Sj)  (COS0Q  -SiSin0Q) 

(8.5) 

X  1  £ 

i 

!  '1 

2 

=  -(I+S2)  (cos0ij-S2Sin0ij)  '  .  ,  ' 

(8.6) 

The  quantity  P  is 

the  joint  probability  density  associated  with  this 

process. 

We,  now  change  va»iab,les  from  (5^,52)  to  (6jj»02R^’  equations  of  trans¬ 
formation  are  (from  Appendix  C  and  Reference  4)  '  ' 


'  30 


(sine^-nsine^y) 

(cos9Q-ncos8j^j) 


So  =  tan 


(TT-e  IT"  02^)1 


(8.7) 


(8.8) 


The  Jacobian  of  the  transformation  is 


0(e  6,,)  '  =  n  3e,2p’'-°lT-°2R)' 

■  ^  (cbseQ-ncoseTT)^  t  2 


(8.9) 


A  little  geometry  shows  that  the  allowed  range  of  6,|^  is  -je^^  <  lOjj!* 

Thus,  since  le^jl  ‘'S  always  less  than  ir/Z,  we  have  102^1  <  it/2,  nd  the  step 

function  in  Eq;  (8.1)  is  always  unity.  The  range  of  Bjy  is  determined  by  the  now 
familiar  conditions  (nj*2)  >  0  and  (-ic^'n,)  >  0.  The  result  is  x**  <  ©ii-  <  X^» 
where 


II  1  /cose„ 

•■=  it/2  -  sin-^  L-5. 


X^  =  00  ■  [^/2-sin"'(^)] 


(8.10) 


(8.11) 


By  use  of  the  transformation  equations,  the  final  form  for  can  now  be  written 
down: 

/  U 

fX 

{tjj^(S,0Q)^  =  (sec0Q)jd0j,T  J(9j^j»6)u(  l0ji  l-l0|)tTR(s2(0TT)  >52(0 >  e)] 

,  , 

X  Pts^vOTj) ,S2(9jt»6) »^j(8o^  *^12^0 IT ^  ^ 


X  exp{-K[)l2(6ij)-X'4)sec0TT} 

X  exp[-K(Jl+£-^)sec0]  , 


where 


n  r(TT-9,T"e)l 

S2(eiT.9)  =  tanj^ - - - J 


(8.12) 


(8.13) 
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The  probability  density  P  can  be  written  as  the  product  of  five  probability 
functions: 


=  P{Si)PiP^P^P^ 


(8.14) 


where 


Pi  = 

PtIl(9o)|Sil  . 

(7.16) 

P2  = 

(8.15) 

P'  = 
^3 

P (^2 1 ^ 12^^1T^ ^  , 

(8.16) 

P '  = 

P [ 12(0) 1 S j ,S2 ,Ii(6o^ >" 12(®1T^ ^  . 

(8.17) 

relations 

for  these  probabilities  are 

(i+2Bq) 

(8.18) 

P2  =  1 


^(®0’®1tI^1^ 


(8.19) 


where  the  bistatic  shadowing  probability  is  given  by  Eq.  (6.9); 

2 


P'  = 


1-erf  V 


IT 


1-erf  I 


u(S2-cot0j.|.)P(S2)  .  O<0jj<Tr/2 


u(cotGjy-S2)P(s2)  »  ejj<0  ,  (8.20) 


P4  = 


1+2B 


(8.21) 


The  quantities  B  and  B^  are  given  by  Eqs.  (5.27)  and  (6.2),  while 


cote 


IT' 


'IT 


(8.22) 
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The  remaining  quantities  which  must  be  specified  are  Jig  and  £4,  the  average 
illuminated  heights  of  points  1  and  2  for  this  process.  The  results  are 


£3(0it) 


cj 


r(U2B„)r(2-^2B„) 

[r(2+2r3^+2Bj.j.)  •  r(l+2Bj^)r(2+2Bjj)] 

=  ^4  =  (i5b)  ^  ^  ■ 


(8.23) 

(8.24) 


IX.  RECIPROCITY,  ANU  TRANSMISSION  INTO  A  LESS  REFRACTIVE  MEDIUM 

In  the  previous  sections,  our  investigations  were  concerned  with  th’  trans¬ 
mission  of  light  from  a  medium  of  index  unity  into  a  medium  of  refractive  index 
n  >  1.  We  are,  of  course,  interested  in  the  reverse  problem  of  transmission  into 
the  less  refractive  medium.  Thus,  for  the  sake  of  argument,  we  could  consider  a 
source  located  at  depth  z  =  -£,  with  plane  radiation  emitted  into  direction 
towards  the  surface.  IV  strength  of  the  radiation  refracted  into  the  air  and 
traveling  in  direction  k'  away  from  the  surface  is  then  of  interest.  Therefore, 

we  want  to  derive  {t''^(kik  )}  . 

'  ^  'n->-l 

It  is  expected  that  the  solution  of  the  (n->l)  transmission  problem  will 
follow  from  the  already  derived  (l-»-n)  transmission  coefficients.  This  follows 
because  the  transmission  coefficient  nust  have  a  reciprocity  property  associated 
with  it  which  relates  to  L  .e  inte»’charige  of  source  and  observer.  Reciprocity  is 
a  general  statement  of  the  dynamical  reversibility  of  the  system,  so  its  matli- 
ematical  expression  must  be  independent  of  the  shape  of  the  surface,  Thev'efore, 
we  will  derive  the  reciprocity  theorem  for  a  plane  surface  and  then  simply  infer 
its  form  for  the  general  rough  surface.  We  will  also  assume  n  to  be  real,  which 
is  quite  valid  for  the  optical  transmission  problem. 

Consider  plane-wave  radiation  incident  on  a  flat  surface  with  normal  in  the 
direction  of  the  z  axis.  We  choose  this  orientation  because  the  transmission  and 
scattering  properties  of  any  stationary  random  rough  surface  must  reduce  in  the 

(  ) 
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long-wavelength  limit  to  the  transmission  and  scattering  results  for  this  plane 
surface.  We  will  initially  restrict  the  incident  polarization  to  be  either  hori¬ 
zontal  cr  vertical  relative  to  the  local  plane  of  incidence;  this  will  provide  a 
reciprocity  relation  for  one  dimensional  surfaces.  The  transmission  coefficients 
of  interest  can  be  written  down  by  inspection  for  the  plane  surface: 

tp(k,ko)  =  [l-r(coseQ)]6^0-sin"^(--||-^!j  ,  (9.1) 

tp(-ko,-l<)  =  [l-f(cos0) ]6[0Q-sin"^{nsine)]  ,  (9.2) 


Where  the  reflectivities  are  those  appropriate  to  either  horizontal  or  vertical 
polarization.  For  example,  Eq.  (9.1)  states  that  the  transmission  coefficient 
is  zero  unless  the  angle  of  observation  coincides  with  the  angle  of  refraction 
for  a  plane  surface  as  calculated  by  Snell's  law.  The  bar  on  the  reflectivity 
in  Eq.  (9.2)  is  a  reminder  that  for  the  (1-^n)  process,  we  have,  e.g.. 


R^(cos0q)  = 


but  for  the  (rh>l)  process. 


cosOq  - 


[TJ - K"~~  'p 

n^-s1n^0J  * 


R^(cos0) 


COS0  -  ^^2  ~  sin^oj 

1  .  2 
K  -  sin'^O 

.XT 


COS0  + 


We  now  use  the  relation 


6[f(0)]  = 

|f’(e*)| 


where 


f(0)  =  0p-sin'^(nsin0) 


(9.3) 


(9.4) 


(9.5) 


(9.6) 
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(9.7) 

(9.8) 

(9.9) 

(9.10) 

(9.11) 

(9.12) 


tp(-ko.-k) 


cose.i  H- 

■  ■  -  t  ( k  k  ) 

ncosej 


(9.13) 


This  is  not  the  most  general  form  of  the  reciprocity  relation,  for  the  t's 
are  defined  "per  unit  sector  of  angle,"  and  angular  widths  are  not  invariant  in 
the  transmission  process.  In  order  to  extend  Eq.  (9.13),  we  again  illuminate  the 
flat  plane  with  radiation  of  constant  intensity  but  now  with  a  narrow  sector  de^ 
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of  illumination  angles  about  the  angle  9^.  From  Snell's  law,  sine^  =  nsine,  the 
angular  width  of  the  corresponding  refracted  radiation  is 


de 


cose^ 
ncos'e  ‘^^o 


The  substitution  of  this  result  into  Eq.  (9.13)  yields 


de  tp(k,kQ) 


(9.14) 


(9.15) 


The  left-hand  side  of  Eq.  (9.15)  is  the  fraction  of  the  incident  energy  which  is 
transmitted.  This  is  evidently  equal  to  the  fraction  transmitted  in  the  recipro' 
cal  process,  i.e.,  with  the  source  in  the  medium  and  radiation  of  width  de  about 
direction  -k  incident  at  the  interface. 

The  value  of  Eq.  (9.15)  is  that  we  can  generalize  to  "two  dimensional" 
results  by  simply  replacing  angular  sectors  by  solid  angles  and  by  introducing 
polarization  notation  in  the  usual  way;  i.e.. 


(d(pdesine)tj‘'(k,kQ)  =  (d(|)Qde^sinep)  tJ'^(-kQ,-k)  .  (9.16) 


Since  the  incident  ray,  the  normal,  and  the  refracted  ray  all  lie  in  the  same 
plane,  we  have  d(|)  =  dcl)^.  Then,  by  use  of  Snell's  law  and  Eq.  (9.14), 


(•'k^j,-k) 


cose„ 

2 

n  cose 


l-Ji, 

P 


t^'(k.,k  ) 


(9.17) 


Because  of  the  general  nature  of  the  reciprocity  principle,  we  expect  that  the 
coefficients  for  the  general  random  rough  surface  satisfy 


cose^j 

n  cose 


Finally,  we  replace  -k-»-k^,  -k^-^k'  to  obtain 


cos9 

2 

n  cose. 


(9.18) 


(9.19) 
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We  can  new  substitute  the  results  of  Sections  V,  VII,  and  VIII  into  the  right- 
hand  side  of  Eq.  (9.19)  to  obtain  the  corresponding  partial  transmission  coef¬ 
ficients  for  refraction  from  water  into  air. 

X.  INTENSITY  RELATION  FOR  UPWELLIN6  RADIATION 

Upwelling  radiation  from  beneath  the  surface  is  an  important  component  of 
sea  color.  Such  radiation  originates  from  the  backscattering  of  transmitted  sun¬ 
light  by  both  molecular  and  particulate  matter  or,  in  shallow  water,  by  the  sea 
bottom.  Some  of  this  scattered  radiation  is  reflected  back  into  the  medium  by 
the  sea  surface,  but  a  certain  fraction  will  be  transmitted  into  the  atmosphere. 
The  intensity  reaching  an  observer  above  the  surface  in  direction  k  will  depend 
on  the  original  intensity,  the  transmission  coefficients  for  the  (rough)  air- 
water  and  water-air  interfaces,  and  on  the  scattering  coefficient  defined  as 


Y^^(ic2»^l)  =  fraction  of  incident  (kj)  radiation,  with 
polarization  i,  which  is  scattered  into 
direction  k2,  in  polarization  state  j. 

The  expression  for  the  intensity  of  the  upwelling  radiation,  in  polarization  state 
n,  is  clearly 

^  i,  1  k  k  k 

j=i  m=l  •’ 


X  t'J(k^,k^,)(cos0^j)I^(k^j) 


(10.1) 


where  I^(kQ)  is  the  initial  intensity,  in  polarization  state  i. 

Strictly  speaking,  the  transmission  coefficients  for  the  physical -ocean 
medium  will  differ  from  those  derived  in  the  previous  sections  on  the  basis  of  an 
idealized  semi -infinite  homogeneous  medium.  The  physical  medium  has  a  fluctuating 
index  of  refraction  which  makes  the  description  of  propagation  much  more  compli¬ 
cated  than  the  exponential  absorption  model  adopted  here.  Nevertheless,  the  t^^ 
of  this  report  should  be  sufficiently  accurate  for  clear  waters.  Also,  we  will 
assume  that  the  distance  to  the  bottom  is  large  compared  with  the  absorption  length 
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so  that  any  radiation  scattered  off  the  bottom  and  subsequently  scattered  back 
into  the  medium  by  the  surface  will  be  of  negligible  intensity  as  compared  to  the 
direct  transmission  component. 

The  simplest  application  of  Eq.  (10.1)  is  to  a  homogeneous  (but  absorbing) 
medium  with  a  plane  bottom  (of  zero  slope).  Then,  the  scattering  coefficient  is 

y-^^ks.kt^  =  rj(k^*z)6^'"6(k5-kp  ,  (10.2) 


where 


ks  =  k^-2(k^*z)z  .  (10.3) 

Here,  r.  is  the  reflectivity  of  the  bottom  for  radiation  in  polarization  state  j. 
We  also  assume  a  plane  wave  for  the  incident  radiation,  so 


(10.4) 


Then,  the  substitution  of  Eqs.  (10.2)  and  (10.4)  into  Eq.  (10.1)  yields 

J  • 


All  quantities  in  this  equation  are  known,  with  the  transmission  coefficients 
available  from  Sections  V  and  IX. 


XI.  CONCLUSIONS 

A  ray  theory  which  describes  the  transmission  of  light  across  a  random  rough 
boundary  has  been  formulated.  The  approach  uses  a  complete  geometrical -optics 
treatment  in  that  multiple-scatter  and  shadowing  effects  are  retained.  The  main 
results  of  the  analysis  are  as  follows: 

1)  A  theory  of  shadowing  for  the  transmission  problem  las  been  derived; 
i.e.,  the  probability  that  a  surface  point  is  visible  to  a  source  above  the  sur¬ 
face  as  well  as  to  an  observer  below  the  surface  (or  vice  versa)  is  now  available 
in  closed  form. 
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2)  The  transmission  coefficient,  which  de^-ermines  the  angular  profile  of 
transmitted  energy,  is  expressed  as  an  expansion  in  terms  of  the  number  of  ray 
intersections  with  the  surface.  All  terms  in  the  expansion  are  positive,  so 
the  retention  of,  ''.ay,  the  first  m  terms  yields  a  lower  bound  to  the  true  value 
of  the  transmission  coefficient.  The  leading  term,  the  direct  transmission  coef¬ 
ficient,  is  corrected  for  shadowing  but  contains  no  multiple-scatter  infoiination. 
As  expected,  it  is  proportional  to  the  probability  density  of  that  slope  which 
connects  the  incident  and  observer  directions  by  Snell's  law.  The  direct  trans¬ 
mission  coefficient  is  derived  for  the  general  two  dimensional  random  rough  sur¬ 
face,  and  it  therefore  contains  polarization  information.  Expressions  for  both 
air-to-water  and  water-to-air  transmission  are  available,  the  latter  by  a  reci¬ 
procity  argument. 

3)  Higher-order  transmission  terns  were  formulated  only  for  the  mathemati¬ 
cally  tractable  one  dimensional  rough  surface  (cylindrical  symmetry).  Thus,  a 
ray  of  sunlight  may  reflect  from  one  point  on  the  ocean  surface  only  to  intersect 
the  surface  again  at  a  second  pjint.  The  angular  profile  due  to  transmission  of 
energy  at  the  second  point  has  been  derived;  the  total  fraction  of  incident  energy 
associated  with  this  double -scatter  process  is  about  3%  of  the  directly  trans¬ 
mitted  energy  for  likely  ocean  roughnesses.  The  corresponding  process  in  which 

a  transmitted  ray  of  sunlight  is  reflected  by  the  under-sea  surface  was  also 
formulated.  K  is  expected  to  contribute  significantly  to  vhe  below-surface 
radiation  profile  only  for  near-grazing  incidence  and  the  very  roughest  of  ocean 
surfaces. 


U 
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APPENDIX  A 


% 

We  consider  the  general  case  of  a  medium  with  complex  index  of  refraction  n, 


where 


'h 

n 


n  +  in- 


(A.l) 


Snell '  nas  the  familiar  fonn  for  radiation  incident  on  a  plane  surface  at 
angle  0^, 

^  sine. 

sine,  =  -r —  .  (A. 2) 

^  n 


but  since  n  is  complex,  the  quantity  Op  is  also  complex  and  no  longe»  has  the 
simple  significance  of  the  angle  of  refraction.  The  angle  that  is  significant, 
e.r,  relates  the  direction  of  propagation  of  the  surfaces  of  constant  real  phase 


1 

with  the  normal 

to  the  boundary. 

O 

From  Born  and  Wolf,  we  have 

sine 

sirej  g 

(A.3) 

where 

D  = 

f  9  7  9 

[Sin  e^+nVC COSY  -  —  sinY)  j 

(A.4) 

and  q  and  -y  are  real  numbers  defined  by 

% 

cos0y  =  q  exp(iY) 


For  optical  radiation  in  the  blue-green  band,  the  actenuation  length  for 

-4  “1 

clear  water  ranges  up  to  about  30  meters,  so  at  the  peak  K  =  kn-  =  3x10  cm  , 

-9  * 

and  is  a  few  tiaies  10  •  When  n^/n  «  1,  the  quantity  D  reduces  to 

D  n|l  -  ^  sin^e^j  j  .  (A. 5) 

-18 

Since  the  correction  term  is  only  of  the  order  of  10  for  the  optical  (blue- 
green)  transmission  problem  (n  =  1.33),  Eq.  (A.3i  takes  the  form 

sine 

sine-  =  ■  ■  ^  .  (A. 6) 

i  n 
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APPENDIX  B 
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with 


Cj  = 


(B.15) 


Cg  -  hj*Py 


(B.16) 


Equa^tion  (B.8)  now  takes  tl;e  form 


=  (ajCjt,  +  bjC2t„)  p,,  e 


ikjT-r 


•1^ 

ikj^.r 


(B.17) 


The  time  averaged  Poynting  vector,  N  =  |-  Re{E  x  H*),  can  now  be  used  to 
find  the  fraction  of  energy  incident  on  the  plane  surface,  in  polarization  state 
H,  which  is  refracted  into  final  polarization  state  H  or  V.  We  designate  these 

uu  uw 

fractions  as,  t  and  t  : 


.  I^frl  /-“I 


•^inc'  ^ 


2|cosaj 


and,  similarly. 


-  njbjCjt^  -  ajC2ty| 


HV  PI 

T  =  n|ajCjty  +  bjC2 

\coscxi, 


where,  by  Snell's  law. 


9 


ncosa 


j  =  [{n^-D  +  cos'^aj] 


(B.18) 


(B.19) 


(B.20) 


VH  VV 

The  alternate  quantities  t  and  t  ,  corresponding  to  vertically  polarized  in¬ 
cident  radiation,  follow  from  Eqs.  (8.19)  and  (B.20)  by  the  replacement 

bj^-^-aj*. 


cosa 


APPENDIX  C 


Snell's  law  relates  the  angle  of  refraction  to  the  angle  of  incidence  on  a 
plane  surface,  but  the  angles  are  defined  relative  to  the  surface  normal.  For 
an  irregular  boundary,  we  must  apply  Snell's  law  to  each  of  the  locally  flat  and 
arbitrarily  oriented  elements  comprising  the  surface.  Then  it  is  convenient  to 
talk  about  the  angles  of  incidence  and  refraction  defined  relative  to  a  single 
coordinate  system. 


For  isotropic  media,  the  incident  ray,  the  refracted  ray,  and  the  normal  to 
thi  surface  all  lie  in  the  same  plane  (Fig.  6).  We  want  to  derive  the  vector 
relation  between  the  direction  of  the  refracted  ray  and  the  directions  of  the 
incident  ):;ay  and  surface  normal.  One  triplet  of  unit  vectors  is 


(kpX 


» 


y\ 

,  and  V,  where 


V 


(V  hi) 

IkpX  kjjl 


(C.l) 


Then,  the  application  of  vector  algebra  leads  to 

A  A  A 

nj^  =  cosa(-kg)  +  sina  V  (C.2) 


By  use  of  Snell's  law,  sino.  =  nsina',  Eq.  {C.3)  simplifies  to 

A  A  A  O  ^ 

nk^T  "  a]  -  cosa}n^ 


(C.3) 


(C.4) 


We  now  apply  the  general  vector  relation,  Eq.  (C.4),  to  the  coordinate 

s 

system  of  Fig.  5.  The  components  of  k^  and  k^j  are 

A  A  A 

kg  =  "sine^^  X  -  cose^  z  ,  (C.5) 
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*^1T  “  ~  (sin9jyCos<!>2j  x  +  sine,jsin(|>jj  y  +  cose^y  z) 


(C.6) 


where  the  minus  sign  in  Eq,  (C.6)  arises  because  kjj  is  defined  relative  to  the  -z 
direction  (Fig.  2).  In  addition,  we  make  use  of  the  geometrical  relation  between 
the  direction  of  the  normal  to  a  surface  point  and  the  slope  components  at  that 
point: 


o  p  ^ 


(C.7) 


The  substitution  of  Eqs.  (C.5)-(C.7)  into  Eq.  (C.4)  yields  expressions  for 
(SjxjSjy)  in  terms  of  the  incident  and  refracted  angles: 


’lx 


’ly 


- 0  - 11 - ^11' 

> 

(C.8) 

(ncose^j-cose^) 

nsinejjSin(l)2j 

(C.9) 

(ncose^j-cose^) 

("'f 

In  Section  V,  it  is  necessary  to  transform  integrations  over  slope  var— 

ables  into  integrals  over  the  angles  (02-j-,(j)jj).  The  corresponding  limits  of 

integration  for  the  angle  variables  are  not  simply  arrived  at  for  the  general 

rough  surface.  We  see,  from  Eq.  (C.4),  that,  for  fixed  k.  and  n,  the  direction 

^  0 

kjy  varies  with  the  normal  n^.  Hence  restrictions  on  the  surface  normal  are 
equivalent  to  restrictions  on  “it-  The  first  restriction  is  n^'Z  >  0.  From  Eq. 
(C.4)-(C.6),  we  see  that  this  condition  is  equivalent  to 


(ncose^j  -  cos9j^) 


} 


>  0 


(n^-l)+cos^a]  -co‘aj 
or,  because  the  denominator  is  always  greater  than  zero. 


ncose 


IT 


cose. 


(C.IO) 
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The  second  condition  is  (-kg'n^)  >  0.  The  substitution  of  Eqs.  (C.5)  and  (C.7) 
into  this  condition  yields 

Six  <  cote^  .  (C.ll) 

We  now  substitute  Eq.  (C.8)  into  the  left-hand  side  of  the  inequality.  After 
algebra,  the  condition  reduces  to 


Since  n 


=  3/4  for  water,  this  latter  condition  is  a  severe  one. 


{C.12) 


i. 
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